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We consider the problem of a single particle interacting with A*' identical fermions, at zero tem- 
perature and in one dimension. We calculate the binding energy as well as the effective mass of the 
single particle. We use an approximate method developed in the three-dimensional case, where the 
Hilbert space for the excited states of the A'^ fermions is restricted to have at most two particle- 
hole pairs. When the mass of the single particle is equal to the fermion mass, we find an excellent 
agreement with the exact results of McGuire. When the mass of the single particle is infinite, we 
solve exactly the problem and find again excellent agreement between approximate results and exact 
ones. This overall agreement in one dimension gives a strong validation for the approximate method 
applied in three dimensions. Moreover it shows that our approximate treatment is excellent for the 
one dimensional problem in the general case with respect to the mass of the single particle. 
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I. INTRODUCTION 

The field of ultracold fermionic atoms has seen during 
the last ten years an impressive amount of experimen- 
tal and theoretical activity These systems provide 
experimental realization of various many-body systems 
that for a long time have been considered as purely the- 
oretical models 0, For instance, the ID Fermi gas 
can be experimentally realized by tightly confining the 
atomic cloud in the radial directions and weakly con- 
fining it along the axial direction. Its behavior can be 
characterized to a very good approximation by an ef- 
fective ID coupling constant g [j]. One major experi- 
mental achievement of the last few years was the realiza- 
tion of the BEC-BCS crossover by driving a mixture of 
two fermionic populations with an attractive interaction 
through a Feshbach resonance. On the BEG side of the 
crossover, a bound state appears between particles of dif- 
ferent species, while on the BGS side, pairs exist only in 
the presence of the whole Fermi sea. In one dimension, 
the situation is quite different, because any attractive in- 
teraction produces a bound state. However, an analog 
of the 3D crossover can be achieved in a quasi-lD trap 
where one can vary the effective coupling constant g from 
—00 to -|-oo by detuning the 3D scattering length 

Recently attention has focused on imbalanced systems 
and in particular on the strongly imbalanced regime. The 
extreme situation corresponds to the case of a single |- 
spin particle interacting with a "f-spin Fermi sea. The ef- 
fective quasiparticle parameters, namely, the binding en- 
ergy Eb of the |-spin atom and its effective mass m* , are 
known to be of major interest for the analysis of the phase 
diagram [1, 0|- In the three-dimensional case, analytical 
and MG calculations have been done [1, 0, i, i, E^JilL 
Quite surprisingly a simple variational calculation Q 
gives results remarkably close to MG calculations. Very 
recently we have explained this puzzling feature and 
shown that this simple calculation is just the first step in 
a series of approximations converging extremely rapidly 



toward the exact result. These approximations corre- 
spond to restrict the Hilbert space to states with in- 
creasing number of particle-hole excitations. The very 
fast convergence is due to almost perfect destructive in- 
terferences occurring when more than one particle-hole 
excitation is present. As a result [T^ a variational calcu- 
lation taking into account at most two particle-hole exci- 
tations provides already an essentially exact solution to 
the problem. In this paper, we use the same approxima- 
tion to address the ID problem, having mostly in mind 
the case of attractive interaction. However wc consider 
also the case of repulsive interaction. In the case of equal 
masses, we can compare our results with the exact solu- 
tion obtained by McGuire 01- The excellent agreement 
we find provides a further check of the reliability of our 
approach. 

II. GENERAL CALCULATION 

Wc consider the problem of a single ^-spin particle with 
mass in the presence of a ID Fermi sea, with Fermi 
momentum kp and density n| — kp/i:, of "f-spins with 
mass TO|. In the limit of short-range interactions, we can 
write the Hamiltonian of the system as: 

H = Hc + V (1) 
iJe = ^£;(Q)6t^6Q+^efccIcfc (2) 

Q k 

V = g ^KK'QQ'C^xCK'h^Q^Q' (3) 

KK'QQ' 

where = A:^/2m|, E{Q) = Q"^ /2mi, and and c\ 
are annihilation and creation operators for "f-spin atoms 
while bk and b\. are for the j-spin atom. The Kronecker 
symbol Skk'qq' ensures momentum conservation in the 
scattering. Gontrary to the three-dimensional case p^ . 
the representation of the interaction potential by a Dirac 
distribution is mathematically well-defined, since there 
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are no ultraviolet divergences. The coupling constant g 
is either negative or positive, and there is no need to 
let it go to zero. It is related to the value of the one- 
dimensional scattering length hy g = — 1/(771^0), where 
nir = mj^m|/(m| is the reduced mass (we have set 

h = 1). 

Following the same idea as in Ref. [l^, we write, for a 
system of total momentum p, the many-body wave func- 
tion as: 

1^) = aobl\0) + Y,akqbl^^_Acq\0) (4) 

kq 



+ 1 H (^kk'qq'bl 



p+q+q' -k-k'^k^k'^q^l' 



44,C,C,, |0) 



kk' qq' 



where |0) = Il|5|<fcF ^9 is the noninteracting Fermi 
sea of t-spins and the sums on q and k are implicitly 
limited to j^l < kp and |fc| > kp- This state is a super- 
position of states with up to two particle- hole pairs. In 
the first term, the |-spin free Fermi sea is in its ground 
state and the ^-spin atom carries the momentum p, while 
the two other terms correspond to excited states with 
creation of one or two particle-hole pairs in the Fermi 
sea, the i-spin atom carrying the rest of the momentum. 
In the third term, the coefficients akk'qq' are antisym- 
metric with respect to the exchange of their arguments 

{akk'qqi = -Ctk'kqq' = -Oikk'q'q = Oik' kq' q) and thc factOr 

1 /4 corrects multiple counting. 

Now, we write the Schrodinger equation H\iIj) = i?|V') 
and project it on the subspacc corresponding to Eq.(|3]). 
This yields a set of three coupled equations. Thc two 
first equations are obtained by projecting onto thc full 
Fermi sea and the Fermi sea with a single particle-hole 
pair: 



-g ^E'^^'ao^^a/cg 



(5) 



kq 



-9 ^E^p]lg<^kq^ao+^aKq -^OkQ -^akKqQ (6) 

Q KQ 



K 



where we have introduced E, 

eH = -E 

p,kq 



(0) 



-E 



p — -1^ -r + E(j)), and 
.grij + Ei(p + q - k) + ek - eg. Here E 
denotes the energy measured from the energy of the free 
Fermi sea. Note that, if we omit the last term of Eq.®, 
these two first equations would be decoupled from the 
third one and this would correspond to the variational 
calculation Q with a single particle-hole pair. The next 
equation is: 

1 (2) 

E^p lf^,^^,akk'qq'^~akq- ak'q'+ akq'+ ak'q (7) 

aKk'qq' +^ O-kKqq' ^ ^ O-kk'Qq' " ^ OLkk' qQ 
K K Q Q 



where E, 



(2) 



p^kk' qq' 



-E + gm;+E{p + q + q' -k-k') + ek + 
Here, since we have restricted the Hilbert 



omitted thc overlap with states containing three particle- 
hole pairs. This restriction corresponds to perform a vari- 
ational calculation. For p = 0, we have E = —Eh, while 
the variation of E for small p gives the effective mass. 
We consider first the calculation of the binding energy. 

It is first worth noticing that, for attractive interaction, 
the two limiting cases, namely g 0- and g —00, are 
already well described by the first order approximation, 
corresponding to retain only Eq.([5]) and Eq.®, with the 
last term in the right-hand side of Eq.® omitted. In the 
weak coupling limit g Eq.Q gives to lowest or- 

der the expected result, that is the mean- field interaction 
energy E = gn-\. Then Eq.® provides the correction 
expected from second order perturbation theory, namely 
E = .gn| - m|.g2[l/8 + (^2(0;) - Li2(-a))/(27r2)], where 
a = (to| — m|)/(mx + m-^) and Li2(z) = J2n>i^"'/''^'^ 
is the dilogarithm function. In thc case of equal masses, 
the result is identical to the one obtained by McGuire 
in the same approximation. 

In the strong coupling regime, we can check that the 
binding energy of the molecular state —E = i^t = 
mrg^/2 is indeed solution to dominant order. In this 
regime, since q is limited to |g| < kp, thc dependence 
on q is negligible in E^^j.^. Now we can sum Eq.® 
over q. This gives -g^^ [Eb + k'^ /{2mr)]J2q (^kq = 
a^kp /T^+^Kq o^Kq- Then we can divide this equation by 
\Eb + A:^/(2mr)] and sum over k. To dominant order we 
find Y.k [Eb + k'^ / [Imr)]^^ = ^Jmr/{2Eb), leading to 

1 + g-^^2Eb/mr J^kq ^kq 

g-^yj2Eb/mr 



-aokp/n. Inserting this 

relation in Eq.® gives 1 -I- g~^ ^/2EiJrri^ = gkp/{TrEb) 
to dominant order. This is indeed consistent with Eb ~ 
mrg^/2. 

In the general case, these equations can be numerically 
solved by successive iterations to provide the binding en- 
ergy Eb and the effective mass m*. This is done in the 
same way as for the three dimensional treatment, but 
obviously the numerical work is much lighter since one 
needs to perform only one-dimensional integrals. On the 
other hand, in contrast with thc 3D situation, the third 
term in the right-hand side of Eq.® and the last two 
terms in the right-hand side of Eq. ([7]) do not disappear. 
However since we use an iterative method this does not 
make a major problem. 



III. EQUAL MASSES 

We first investigate the situation of equal masses m| = 
mi = m. In this case, the exact solution of this prob- 
lem has been obtained by McGuire using the Bethe 
ansatz (with a small sign mistake in the last two equa- 
tions of page 126 and in the expression for m*). With 
our notations (we have the opposite sign convention for 
the interaction) thc correct expressions are: 



space to contain at most two particle-hole pairs, we have 



Ep 



-y^ + (1 + y^) arctany 



(8) 
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FIG. 1: (Color online) Reduced j-atom binding energy Eb as 
a function of mg/2kF for equal masses = m| = m and 
attractive interactions g < 0. From top to bottom: numeri- 
cal results for the first (dashed-dotted line) and second (solid 
thick line) order approximations and exact solution Eq.(|8]) 
(solid red line) . The inset shows a zoom of the region around 
mg/2kF = — 2 (we have chosen a value larger than the one 
discussed in the text to make it easier to see the difference on 
the figure). 



m 
m 



[l — - arctany) 



1-2 



arctan y 



l+y' 



(9) 



where y = mg/2kp and Ep = fc|,/2m is the Fermi energy 
of the t atoms. They allow us to check the quality of 
the convergence of our successive approximations to the 
exact many-body solution. 

We begin by considering the case of attractive inter- 
actions g < 0, which is the more interesting one in re- 
lation with the 3D case (T^]. The overall result of the 
comparison is displayed on Fig[TJ As it was indicated 
in Ref.[§] and is clearly seen in Fig[Tl the lowest order 
approximation for the binding energy is already in good 
agreement with the exact solution and the two curves are 
barely distinguishable on the graph. However, the sec- 
ond order approximation improves it significantly, lead- 
ing to results completely undistinguishable from the ex- 
act result on the graph. To be more specific, let us 
consider the case for mg/kp = —1- The exact result 
is El, = 0.9373 E'f- We may consider the Hartrece ap- 
proximation Eh = —gn-^ = —{2/T:)(mg/kp) Ep as a ze- 
roth order approximation. In the present case it gives 
Eb = 0.6366 E'f, which means a 3.2 10~^ relative error 
compared to the exact result. The first order approxi- 
mation gives Eb = 0.9242 Ep, that is a 1.4 10"^ relative 
error. Finally we find for our second order approximation 
Eb = 0.9361 E'f, which corresponds to only a 1.3 lO"'^ 
relative error. We see that the convergence is very fast, 
which gives very strong support to the idea that the rc- 



FIG. 2: (Color online) Reduced J,-atom binding energy Et, as 
a function of mg/2kF for equal masses m| = = m. Same 
conventions as in Fig[l]are used. 



suit obtained in the three-dimensional case is nearly the 
exact solution. Actually this is even more so because the 
convergence in ID is expected to be slower than in 3D. In- 
deed the rapidity of this convergence is related [l^ to the 
effective phase space corresponding to particle wavevec- 
tors k > kp, compared to the phase space corresponding 
to the hole wavevectors q < kp. The ratio between the 
volumes of these phase spaces is clearly an increasing 
function of space dimensionality. 

Naturally, for large values of the coupling constant, all 
the results look very close on FiglT] because they are near 
the large asymptotic value —eb/Ep = — 2 {rng /2kpY . 
In order to make the differences easier to see, we have 
plotted in Fig[2] the difference between the results and 
this asymptotic expression. 

We may also apply the same cascade of successive 
approximations as Ref . (l^ . This leads to the approxi- 



mations E. 



(1) 



E, 



'^Iq in the first order equations and 
in the second 



^p.kk'qq' — ^p.kk'OO ^p,kk'qq' ~ P,kk'q0 

order equations. The results are displayed in FigO for 
mg/kp = —1. 

Finally, the relative effective mass m* /m is plotted 
in FigHl Naturally, in the weak coupling limit, m* ap- 
proaches m since the quasiparticle has the effective mass 
of the free J,-spin atom. In the strong coupling regime, 
TO* approaches the mass 2to of the two bound particles 
forming a dimer. as it is expected physically. The main 
quantitative feature of our results is that, contrary to 
the binding energy, the effective mass is more sensitive 
to our approximations as one could expect, especially in 
the strong coupling regime. The first order approxima- 
tion is substantially improved by our inclusion of two 
particle-hole excitations. As above, this second order ap- 
proximation is fairly close to the exact solution. We note 
that although our variational calculation gives naturally 
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FIG. 3; Reduced J,-atom binding energy Et and efTective mass 
m* , for equal masses m| = mj. = m and mg/kp = — 1. Suc- 
cessive approximations to the exact many-body problem dis- 
cussed in the text are applied. Diamond: first order approx- 
imation. Circle: second order with qq' = 00. Plus: second 
order with qq' = gO. Square: second order with no q approx- 
imation. Cross: exact solution. 




FIG. 4: (Color online) Relative effective mass m* /m as a 
function of mg/2kF for equal masses = = m. Same 
conventions as in Fig[l]are used. 



an upper bound for the energy i?, the obtained effective 
mass is not a lower bound of the exact one since it is 
given by the variation of E. 

It is also of interest to investigate the case of repulsive 
interactions 5 > 0, since it is no more difficult to calcu- 
late than the attractive one. The result for the energy 
is displayed in FigO One sees that for small couplings 
the first and second order approximations are still excel- 
lent, which can merely be understood by continuity from 
the attractive case with small coupling, or the fact that 
in this domain standard perturbation theory is expected 
to be good. However for larger coupling the first order 
approximation becomes markedly less satisfactory. The 
second order approximation remains quite good, but also 




FIG. 5: (Color online) Reduced J,-atom binding energy Et as 
a function of mg/2kp for equal masses m-f = — m and 
repulsive interactions <? > 0. Same conventions as in Fig[T]are 
used. In the strong repulsive limit g —> 00, first and second 
order approximations give E = —Eb = 1.389 i?F (indicated 
by the dashed line) while the exact result is -E = —E^ = Ep. 



with a trend toward a deterioration of the agreement. 
This feature is actually easy to understand. For strongly 
attractive interaction, the large negative energy E makes 
^'pkq ^'^'^ ^'pkk'qq' l^rge and positive. As a result they 
depend weakly on the q (hole) variables, which is the ba- 
sic reason [l^l for the success of the approximations. On 
the other hand for repulsive interactions, E becomes pos- 



itive and E^l^ and i?p7-fc'qq' ^'^'^ smaller, and accordingly 
their relative change under variation of the q variables 
becomes more important. As expected under these cir- 
cumstances, the agreement for the effective mass, shown 
in Figini is very rapidly unsatisfactory for the first order 
approximation. The second order approximation stays 
quite good until mg/2kp ^ 1.5, but the disagreement 
increases quite markedly when the coupling is increased 
much further. These results are interesting since they 
come as a confirmation of the basic reason (l^ for the 
quality of our approximation in the attractive domain. 

In the strong repulsive limit g 00, the binding ener- 
gies given by the first and second order approximations 
are seen to converge to the same value. This limiting 
value can be determined analytically by performing a 1/g 
expansion in the following way (we make the calculation 
directly for a general value of the mass ratio r = mi/m-\ 
since this does not make any problem). We consider the 
equations of the first order approximation, that is Eq.fO 
and Eq. ^ with the last term in the right-hand side omit- 
ted. We set for simplicity ckq = 1 and expand akq in 
powers of 1/y where y = nirg/kp: 



p(2) 



akq 



a 



(10) 



5 



To lowest order, Eq.® and Eq.® gh 

E(o) 



kq 



kq 



(0) 1 



(0) 



(11) 

(12) 



K 



Since in Eq. p^ the sum on K in the right-hand side 
converges, this imphes that we have a 



oo. Taking the hmit k 



(0) 
kq 



for k 



oo in this equation, we obtain 
— 1 and accordingly, from Eg. p^ . a^^^ = 



a), does not depend on q. Note that, after summation 
over q, this finding is completely consistent with Eq. pTj) . 
Now to first order in l/y, we find: 



kp 

ITlr 



(Eb + E{q -k) + ek 

kp 
rrir 



kq 
\ (0) 



(1) 
kq 



-E"S 

K 



E«S 



(1) 



(13) 



(14) 



JO) 



Now we can sum this last equation Eq. (fT4|) over q, 

which gives n-f [Eb + k'^/{2mr) + (1/r - l)Ep/'i) 
Eb, taking into account Eg. p^ . Using the relation 



J") 



-1, we get a closed equation for Eb which 
can be easily solved numerically. For equal masses r = 1, 
it gives E = -Eb 1.389 while from Eq.® the ex- 
act result is i? = —Eb = Ep. One can see that working 
to the second order approximation does not bring any 
change to this result, as supported by numerical results. 
However we will not give details on this since this matter 
is not easy and not a major interest. 




FIG. 6: (Color online) Relative effective mass m*/m as a 
function of Tng/2kF for equal masses m-f — mi = m and 
repulsive interactions g > 0. Same conventions as in Fig[T] 
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IV. UNEQUAL MASSES 

In the following, we come back to the general case 
of unequal masses. We restrict ourselves to the attrac- 
tive case and consider as an example the single value 
mig/kp = — 1 for the coupling constant. The results 
of our calculations as a function of the mass ratio r — 
mi /mi are presented in Fig[7]and[8l respectively for the 
binding energy and the effective mass. For this value of 
the coupling constant the binding energy does not display 
very large variations. The limiting cases can be under- 
stood fairly easily. 

For r — > 0, at fixed coupling constant g, the binding 
energy Eb goes to —gn-^- Indeed in this case, the kinetic 
energy of the |-spin particle goes to infinity as soon as 
its momentum is not zero. This implies that E^^j,^ and 

(2) 

^p.kk'qq' go *o infinity, so that akq and akk'qq' go to zero 
in Eq.® and Eq.([7]). Clearly this argument could be 
extended to any order in our successive approximation 
scheme. Hence this result is exact. This is quite reason- 
able physically since this limiting situation can also be 



0.6 — --L^ — — — — — ^ 
0.001 0.01 0.1 1 10 100 1000 

FIG. 7: (Color online) Reduced j-atom binding energy Eb 
as a function of r = mi/irif for m-\g/kF = —1. Dashed- 
dotted line; first order approximation. Solid thick line: sec- 
ond order approximation. For r ^ oo, Eb goes to (1 + 
1/iv)Ef ~ 1.318 (Eq.^), while for small r, Et/Ep ^ 
2/7r — 2r logr/TT^ (dashed red line). 

seen as the case of a single particle of mass m^ moving in 
the presence of particles of infinite mass, which means 
that there is no dynamics associated with them. Hence 
it is quite natural to find only the mean-field term —gn-\ 
in the binding energy. 

The first correction is obtained from Eq.®. To dom- 
inant order, we find a^q = —ga^j \E[q — k) + ek — e^]. 
Inserting this relation in Eq.(I5]) gives Eb = —gni — 
g^m^rlog r/n'^ . 

As seen from FigjSl the variations of the effective mass 
are very small, which is the essential result. It is not clear 
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FIG. 8: Relative effective mass m* /mi as a function oi r = 
my/m-\ for mig/kp ~ —1. Dashed-dotted line: first order 
approximation. Solid thick line: second order approximation. 
For r and r oo, m* goes to mi. 



how the detailed structure would survive in improved ap- 
proximations. The limiting cases arc again easy to ob- 
tain. For the same reason as above, the effective mass 
in the limit r ^ is merely equal to the bare mass 
m* = mi- And the same is obviously true in the op- 
posite case where r — > oo, which we consider now with 
respect to the binding energy as a function of the cou- 
pling constant. 



I 



-10 



-15 



-20 



-25 



-4.2 



-4.4 



-4.6 




-2.05 



-1.95: 



-30 

-5 -4 -3 -2 -1 

FIG. 9: (Color online) Reduced j-atom binding energy Et 
as a function of m-ig/kp for mi = oo and attractive interac- 
tions g < 0. From top to bottom: numerical results for the 
first (dashed-dotted line) and second (solid thick line) order 
approximations and exact solution (Ea. (|19|) ') (solid red line). 
The inset shows a zoom of the region around m-ig/kp = —2. 
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V. INFINITE MASS 

We consider now in more details the other limit where 
the mass of the J, particle goes to infinity —foo. In 
this case, there is no dynamics associated with this par- 
ticle and the problem reduces to the one-body problem 
of an impurity interacting with a free Fermi sea through 
a contact potential. It can be solved exactly in a way 
analogous to the 3D case. The impurity is at the origin 
and the system is enclosed in the segment [—R,R], with 
i? — > oo in the thermodynamic limit. Without loss of 
generality we assume for simplicity that the number of 
t-spin atoms N is even. The only scattered states cor- 
respond to even wave functions (p'^{x), while odd wave 
functions f^{x) do not feel the contact potential and 
behave just as free particles: 

^+{x) = AkCos{k\x\+S{k)) (15) 
ip;:{x) = ^/T/Rsm{kx) (16) 

where the phase shift is given by cot 5{k) — ka and A^. ~ 
1/R to order which is sufficient for our purposes. 

Since the wave function has to be zero for x ~ R, the 
allowed wave vectors k^ and k~ are given by fcj^i? + 



FIG. 10: Reduced J,-atom binding energy Eb for mi = oo and 
mig/kp = —1. Successive approximations to the exact many- 
body problem discussed in the text are applied. Diamond: 
first order approximation. Circle: second order with qq' — 00. 
Plus: second order with qq' — qO. Square: second order with 
no q approximation. Cross: exact solution. 



S{k^) = (p — l/2)-K and k^R = pi:. For positive cou- 
pling constant g, the N "f-spin atoms occupy TV scattering 
states {N/2 odd states and N/2 even states) correspond- 
ing to p = 1, • • • , N/2. Thus, the Fermi momentum is 
given by kpR ~ {N/2)-k and the energy of all the atoms 

is = Ejf? (fc;)V2mT + E^i? (V)'/2"^T- For nega- 
five coupling constant g, our (5-potential also supports 
a bound state with binding energy — l/2m^a? — 
m^g^ /2 and its wave function is 'i)[x) = e~^^^l°' / ^fa. 
Starting from the non-interacting situation, we conclude 
that the N |-spin atoms occupy N /2 — 1 even scatter- 
ing states corresponding to p = 2, • • • , N/2, the single 
bound state and N /2 odd states. Therefore, the energy 

of all the atoms is £; = -\/2m^Q^ + Y!^^!X [k+f /2m^ + 
(V)V2mT. 

We then identify the difference between the non- 
interacting Fermi sea energy and this energy with the 
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FIG. 11: (Color online) Reduced |-atom binding energy Eb 
as a function of m^g/kF for m| — oo. Same conventions as 
in FigU 



|-atoni binding energy. Since the odd states are not scat- 
tered, we find to order 1/A^ : 



Ef 



N/2 



(17) 



where y = niig/kF and 9 is the usual Heaviside func- 
tion. To order 1/A^, the phase shift writes (5(fc+) = 
— cot~^{2p/{yN)). Wc convert now this sum to an in- 
tegral and find: 



y'o{-y) 



u cot 



du. 



Finally, wc obtain: 



(1 + J/^) arctany 



(18) 



(19) 



This result is valid both for positive and negative which 
is not obvious since there is a bound state in the latter 
situation but not in the former one. This situation is the 
same as the one found for equal masses 0- Note also 
that this result is just one half of the exact result for 
equal masses Eq.([8]), provided we take in the definition 
of y the reduced mass, i.e. y = nirg/kp. 

Here again, we can check the convergence of our theo- 
retical scheme [l^. This is done in Fig [51 which is very 
similar to FiglT] Similarly the inset shows a zoom around 
m^g/kp = —2, where it is clearly seen again that, al- 
though the first order result is already quite good, it is 
markedly improved by the second order one. Specifically 
for m-^g/kF = —1, the first and second order approxima- 
tions give respectively Ei, = 1.302 E'j? and Ejj = 1.317 
while the exact solution is Ei,/Ef = 1 + I/tt ~ 1.318. 
We also show in FiglTOl together with these values, the 



FIG. 12: (Color online) Reduced J,-atom binding energy Eb as 
a function oi m^g/kp for = oo and repulsive interactions 
g > Q- In the strong repulsive limit g — > oo, first and second 
order approximations give E = —Eb — 0.5457 Ef while the 
exact result is i5 = ~Eb — 0.5 Ef- Same conventions as in 
Fig|5]are used. 



results for the intermediate approximations, mentioned 
above for equal masses. 

Considering the large coupling constant limit, we re- 
move again from the binding energy the trivial contri- 
bution ef, of the bound state and display the result in 
FigEH As it can be seen from Eq. p^ . the exact re- 
sult is Eb — eb = 0.5 Ef in this limit, with a corrective 
term —{2/3Tr)EF/\y\. Although the relative error is quite 
small in this regime, as it can be seen from Fig[Tll the 
absolute error becomes larger, of order of 0.1, for the 
first as well as for second order approximation. This is 
clearly due to the near singularity present in the equa- 
tions, linked to the existence of the bound state. 

The corresponding results are given in Fig ll2l for pos- 
itive coupling constant. Just as for the case of equal 
masses, the approximations get worst when the coupling 
constant increases, but the deterioration is markedly less 
pronounced. For g — > oo the exact result is easily found 
to he E = —Eb = 0.5 Ef. The variation of the bind- 
ing energy when g goes from -|-oo to — oo is therefore 
eb + Ef. This is easy to understand. The occupied plane 
waves have exactly the same wave vectors in the two lim- 
its (the phase shift is 7r/2 (modTr) in both cases). The 
only difference is that, for g — + — oo we have the bound 
state with binding energy eb, and in addition a "f particle 
has been transferred to this bound state from the Fermi 
level, which explains the additional contribution Ep. It 
is worth noticing that this variation is half the one found 
in the case of equal masses. 

Finally it is easy and of interest to calculate the |-spin 
local density from the above considerations. It is given 
by the sum of absolute squares of the wave functions. 
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FIG. 13: Reduced |-spin local density as a function of kpx for 
various coupling constants y = mig/kp- Solid line: y = —0.5, 
thin dashed line: y = 0, dashed line: y = 0.5, dashed-dotted 
line: y — ±oo (in the case y = — oo there is in addition a 5 
function spike at a:: = 0). 



Converting this sum to an integral yields: 



= 1 -7rye2'=^l^|y6'(-2/) 



(20) 



+ 



/ [cos^ cot ^ (u / y)) — cos'^ (ukpx)] 

Jq 



du. 



For y ~ Q one has naturally n-^{x) = n^, while for y — > 
±00 this result gives n^{x) = ri| [1 — s\n{2kpx) / {2kFx)\ 
which vanishes when x = 0. Note that this exact result 
Eg. ((20)) for infinite mi is significantly different from the 
exact result found for equal masses 0. In this latter 
case, the t-spin atoms density is always lower than n| 
for positive g and greater than ri,| for negative g p^ . 
The resulting n^{x)/n-\ is displayed in Fig ll3l for several 
values of y. 



VI. CONCLUSION 



We have considered in this paper the problem of a sin- 
gle particle ("impurity") in the presence of N identical 
fermions, at zero temperature and in the one-dimensional 
case. The only interaction term is between the impurity 
and the N fermions, which have no interaction between 
themselves as suited for the effective interaction of a sin- 
gle species of ultracold fermions. This situation can eas- 
ily be realized experimentally in ultracold Fermi gases. 
This problem can be considered as an extreme case of 
polarized Fermi gases and its solution is relevant for the 
case of strong polarization. We have used an approx- 
imate method developed in the three-dimensional case, 
where the Hilbert space for the excited states of the N 
fermions is restricted to have at most two particle-hole 
pairs. We have compared our approximate results with 
exact ones, which are known for special values of the 
mass ratio. When the mass of the impurity is equal to 
the mass of the fermions, we find an excellent agreement 
for the binding energy as well as the effective mass in 
the attractive case. For the repulsive case, which has 
not been considered in three dimensions, the agreement 
is less satisfactory, but this easily understood. The other 
exactly soluble case is the one where the mass of the im- 
purity is infinite. The problem reduces to the one-body 
problem of a static impurity interacting with a free Fermi 
sea, which we solve exactly. Once again the agreement 
of our approximate results with the exact ones is excel- 
lent. This overall agreement in the one-dimensional case 
brings a strong validation for our method applied in three 
dimensions. Moreover it provides an excellent approxi- 
mation for the one-dimensional problem in the general 
case where the mass of the impurity is not simply related 
to the fermion mass. 

We arc grateful to X. Lcyronas for stimulating dis- 
cussions. The "Laboratoire de Physique Statistiquc" is 
"Laboratoire associe au Centre National de la Recherche 
Scientifique et aux Universites Paris 6 et Paris 7" . 
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